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1. Introduction 



Chains of random matrices describe in the double scaling limit N — > oo two- 
dimensional gravity cjaupled to matter which in turn is represented by rational con- 
formal field theories a&B. A connection of the one-matrix models with the theory 
Korteweg-de Vries equations was first established in □. In the series of papers 
the identification of the matter fields with the critical Ising model, Lee- Yang edge 
singularity, and the tricritical Ising model was achieved by comparison of the critical 
exponents (field, dimensions) . A systematic and constructive approach to arbitrar y 
matrix models LrB was made possible by shifting the interest to representations of 
a Heisenberg algebra in terms of Gelfand-Dikii quasi-differential operators IH This 
is a special aspect of Korteweg-de Vries equation theory. All this progress in matrix 
model research was obtained by application jof the technique of orthogonal polynomi- 
als which was developed a long time ago IlirtLx. More details of the history of matrix 
models can be found in review articles ie^. 

After the completion of the work I3~a interest concentrated on the verification 
of the universality classes [p,q] where p and q are coprime, in different realizations. 
From the one- matrix models the sequence of classes [2/ + 1,2], I > 1 waSpObtained and 
only thesa. TwoJiiatrix and some three-matrix models were studied in H,in appendix 
C of Q, Qcl and tl~t3. But the universality classes explicitly constructed were not 
many : [4,3], [5,4], and [8,3]. The claim that^all [p,q] classes appear in two-matrix 
models was put forward (see the abstract of a) . This assertion is easily verified by 
the construction proposed by us. 

Our perturbative approach is not based on an investigation of the Heisenberg al- 
gebra. On the contrary we construct the representations of this algebra by a pertur- 
bative solution of the Schwinger- Dyson equations (synonymous: recursion relations). 
Generally high orders of perturbation theory are needed and a systematic treatment 
of this expansion is therfore crucial. 

The Schwinger-Dyson equations (section 2) involve two multiplication operators 
B 1 and B 2 and two differentiation operators A x and A 2 . We postulate that B\ and 
B 2 tend to differential operators in the scaling domain 

oo 

B x -> a~ 1 ^ a~ ni Q {n \ Q m = Q (1) 

ra=0 
oo 

B 2 -> a~ /l7 J2 CL- n ^P {n \ P (0) = P (2) 

n=0 

If Li(L 2 ) is the degree of the potential Vi(V 2 ) in the two-matrix model action, then 
necessarily 

l<k<Li (3) 
All {Q^}, {P^} are differential operators of order 

ordQ (n) =ordQ + n = l 2 + n (4) 

ord P {n) = ord P + n = l x + n (5) 
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We shall always assume h > l 2 . 7 is the "string susceptibility exponent" and is 
negative, 

jf = <> 2 - < 6 > 

and N tends to infinity. 

If the ansatz ([!]), @ is satisfied we speak of a "singularity of type (h,l 2 )" in the 
action of the model. We can choose this singularity "maximal" by setting 

h = L 1 ,l 2 = L 2 (7) 

We are convinced that all types of singularities that exist can be found if we let L%, 
L 2 run over all naturals. 

After solution of the Schwinger-Dyson equations the commutator 

[B 1 ,A 1 ] = [B 2 ,B 1 ] = [A 2 ,B 2 ] (8) 

is diagonal. Setting it equal the unit operator gives one nonlinear differential equation. 
The perturbative expansion starts from a linear system of equations for the deviations 
of B\ , B 2 from their limiting difference operator. This linear sys tern has corank l 2 — 2. 
It uses two functions Ui , v_i as basis functions. Use of biorthogonal systems of 
polynomials introduces a gauge degree of freedom that affects the whole perturbation 
theory. The "susceptibility function" u is gauge invariant 

u = ui + V-i + a~ 27 M 1 v_ 1 (9) 

(this form for u is valid in our preferred normalization, see section 3). The nonva- 
nishing corank l 2 — 2 permits the introduction of l 2 — 2 new gauge invariant basis 
functions {wi} 1 ^ 2 . Higher order perturbation theory gives one nonlinear differential 
equation for each one, and these equations are ident ical to those obtained from the 
commutator (8) if we postulate that it be diagonal. 

After we study several examples in detail we can summarize our results as follows. 
There are two series of singularity types : 

(A) l 2 does not divide l\ (remember l\ > l 2 ). Then only theJeading orders in 
(jl]) , (H) are relevant and in the sense of quasi-differential operators 113 



The susceptibility exponent is 



where the universality class is 



)^P = qY 2 (10) 



-2 , x 

7 = (11) 

1 p+q-1 V ^ 



\p,q] = [hM (12) 
(B) l 2 divides /1, but l 2 7^ 2. Then 

P = Q h ^ (13) 

and 

[B 2 , Bt] -> a ~ {h+h+lh [P, Q] + higher order terms (14) 
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where 

P = pW _ Q M ~ m Q (1) Q (m - 1) (M = £) (15) 

Km<M 



We show that 
with 

and correspondingly that 
with 



P = c{h,h)Qf (16) 

c(/i, fc) = ; (17) 

h 

7 = — 7 18 

p + q-1 



[p,q] = [l 1 + l,l 2 ] (19) 

Some comments are in order in this context. First there is no universality class 
for singularities (2m, 2). We will later in (section 5) give an argument why this is so. 
Next a universality class with 

p = nq + 1 , n G IN (20) 

can be realized by a singularity either in (A) or in (B). Finally universality classes 
where p and q have a common divisor but this divisor is not q itself, are contained in 
(A). The symmetric case l\ = l 2 belongs to series (B). 

We study the double scaling limit with iV — ■> 00 and only one coupling tuned to 
a critical value. Nevertheless the complete set of nonlinear differential equations for 
u and {wi} is derived. Only the "integration constants" are ze ro and x is the single 
variable in the equations that is kept. 

In order to make this work self-contained we introduce matrix models and discuss 
their analysis by the orthogonal polynomial method in section 2. Two elementary and 
easily proved propositions which are not new but play a crucial role in the perturbative 
a pproach are given the form of "Lemmas" . In section 3 we develop the perturbative 
analysis of the Dyson-Schwinger equations systematically, the decisive observation 
being that the corank of the linear approximation is l 2 — 2. This method is complet e 
in the sense that any singularity l 2 ) can be treated with it. In the subsequent two 
sections we study examples from the series (A) : (4,3) , (5,4) and (6,4); and from 
the series (B) : (3,3) and (6,3). The two examples (4,3) and (5,4) (where l\ and 
I2 are coprime) serve only the purpose to prove that our method reproduces exactly 
the known results. Moreover we show in section 5 that singularities (2m, 2) do not 
fulfill the commutator equation. Expectation values, in p articular their role in tests 
of gauge invariance, are studied in section 6. 



2. The orthogonal polynomial approach to matrix models 

Matrix models based on iV x N hermitean or unitary or other matrices as dynam- 
ical variables have a history of about twenty years. It has turned out that a class of 
"solvable" matrix models exist which are built on a finite number r of such matrices 

{M (a) } (21) 
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which are bilinearly coupled to a chain . In the case of hermitean matrix models (to 
which we shall restrict our attention) this means an action 



S(M (1) , M (2) , . . . , M (r) ) = Tr |£ V a {M (a) ) - ]T c a M (a) M (a+1) ] (22) 

lct=l a=l J 

where each V a is a polynomial of degree L a (only L a > 3 is of interest) 

V a (x) = £ ^-x fc (23) 
fc=i fc 



The action is stable if and only if 



L a is even, all a 

» 



(24) 



[ g), is positive, all a 
In the stable case the partition function of the model is 

Z = [ ndM (Q V 5 (25) 



dM {a) = n^ReM^n^ImM^) (26) 

i<j k<l 

If the action is symmetric under 

M {a) <-> M (r+1 " a) (27) 

we call the model "^-symmetric" , otherwise "asymmetric". A model may have a 
Z2-symmetric critical point which can be approached in an asymmetric fashion such 
as an Ising model with external field H at the critical point. 

When we claim that such matrix models (§) are "solvable" , we mean solvability in 
the scaling domain of a critical point which implies the existence of a full asymptotic 
series 

oo 

(O) ~ ]T a~^uj m {z) (28) 

m=0 

with a~ 7 —>■ in the scaling domain, z a scale invariant variable and computability of 
all u) m by analytic methods, for any "observable" Q. These observables include ther- 
modynamic objects such as logZ or correlation functions of the "Mehta-accessible" 
class (see section 6). Our aim is to develop an algorithm so far, that for a given action 
all singularities can be extracted, the universality class can be described for each case 
and non-leading orders in the scaling behaviour can be calculated to any precision 
desired. 

The full richness of structure develops only in asymmetric two-matrix and many- 
matrix models, which have not been studied carefully before. The symmetric models 
can often be obtained thereafter by a reduction of degrees of freedom. In any case 
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we apply the method of orthogonal polynomials (by Bessis and Mehta, HjO). By a 
theorem of Mehta the action can be expressed in terms of eigenvalues 

{Aj a) }f of M (a) (29) 



N 

»1 



s({xn)=j: 



i=l 



E^(\ w )-E^ a) Af +1) 

.a=l a=l 



(30) 



The partition function goes into 

Z = C r (N) f f[ d/i(AM)A(A«)A(AM) exp[-S({Al Q) })] (31) 

J a=l 

where the measure is 

N 

dMA W ) = IlA (Q) , (32) 

i=l 

and 

A(A) = II(^ - A i) (33) 

is the Vandermonde determinant. For the solvability it is crucial that both A(A) 
appear in the numerator. 

The biorthogonal polynomials 

{n m (A),n n (/i)}^ n=0 (34) 

are defined by 

deg n m = deg fl m = m (35) 

and 

/ n rfA^n m (A«)fl n (AW)exp{- £ V a (X ia) ) + J2c a X {a) X ia+1) } = S mn (36) 

a=l a=l a=l 

through the Schmidt orthogonalization procedure. But this procedure determines 
only 2m + 1 of the 2m + 2 coefficients in IT m , IT m . In fact we have the 
First Lemma: 
Let 

n m (A) = s m A m + 0(A m - 1 ) (37) 
UM = 5 m /i m + 0(/i m ^) (38) 

then the orthogonalization procedure determines only 

SmSm (39) 

but not s m and s m separately. □ 
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This leads to the following corollary: We can of course require that 

s m = s m > (40) 

as in the Z 2 -symmetric case, but expectation values of observables ought to depend 
only on s m s m and not on such arbitrary definitions. They do, as we shall see, but in 
a nontrivial fashion (related with the gauge invariance of Yang-Mills theory observ- 
ables). 

We follow the usual line and introduce differentiation matrices 

n' m = ^(^!) mn n n (4i) 

n 

fl' m = ^(A r ) nm fl n (42) 

n 

and r multiplication matrices B a 

AIT m = XXfli) ron n„(A) (43) 

n 

liU m = ^2(B r ) nm U n (ii) (44) 

n 

(B 2 , B 3 , ... , B r _i are not given here) 

so that 

[B 1 ,A 1 ] = [A r ,B r ] = l (45) 

By arguments typical for the derivation of Dyson-Schwinger equations one can derive 
a set of "equations of motion" 

A 1 + ClJ B 2 = V'^Bi) 

c a _ x B a _ x + c a B a+1 = V' a (B a ), 2 < a < r - 1 

A~ + c r _i5 r _! = V' r (B r ). (46) 
If these equations are satisfied then 

[B^A-y] = Cl [B 2 ,B x \ 

■ ■ ■ = c a [B a+ i, B a ] — • ■ • 

= ^[BnB,.!] =[A r ,B r ] (47) 

can be shown to hold easily. Moreover the following restrictions on the domain of 
A±, A r , B a , where these matrices have nonzero entries, can be shown to be consistent 
with the Dyson-Schwinger equations: 

r 

(ii) m>n = except for - ]J (L a - 1) < n - m < -1 (48) 

a=l 

r 

{Ar) m , n = except for 1 < n - m < JJ (L a - 1) (49) 

a=l 

r r 

{B a ) m , n = except for - J[ (Lp - 1) < n - m < J] (L p - 1) (50) 
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Our strategy is to solve the Dyson-Schwinger equations (f46|) perturbatively under the 
constraints (|48"D-(|50|) and impose the commutator equations ( |45| ) at the end. This 
approach is considerably simplified by the following 
Second Lemma: 

If ( pEpp and p8|)-(|50|) hold, then the commutators (f|7|) are diagonal. 

To prove this lemma we remind the reader that the chain ( JT7D of equations follows 
from (|46f). On the other hand from (f48|) - fl50f ) one can conclude that 

[Bi, Ai] is a lower left triangular matrix and 

[A r , B r ] is an upper right triangular matrix. 
Equality of both commutators proves that they are diagonal. 

Then we can express the equation (f45[) by 



{B\Ai) nn = Ci(-B 2 -Bi) nn = • • • 

= n + const, indep. of n (51) 

3. The perturbative expansion in the double scaling domain 

Those coupling constants that appear as expansion coefficients in the potentials 
V a (p2[), ( p3[) are fixed to critical values. The corresponding critical action is then 
multiplied with 

* (52) 
9 

and tuned as 

N -> oo, g -> c/ c . (53) 
The matrix labels n, m become continuous variables in this limit 

AT s 

< ^ < 1 (for the labels of the eigenvalues). (54) 
With the help of the string susceptibility exponent 7 (0 < —7 < |) one defines 

i = a 2 ^ (55) 

f = ^£(i- a 2 x) (56) 
9 

g = g c {l-a 2 z) (57) 

z = (g c -g)N^; (58) 

2; is the single scale invariant variable, later we will introduce the "integration con- 
stants" h s that play a similar role, and 7 is as in (|TTf) , (|l8|). 

For the two- matrix models, be they symmetric or not, we proceed as follows: We 

set 

(-Bi) n>n+fc = r k (n) (59) 
(B 2 ) n ,n+k = s k (n) (60) 
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and define 

rk(n - f) = pk + a~ 2l u k {x) (61) 
s k {n-\) = a k + a- 2 ^v k {x). (62) 

The translation of the arguments of r k , s k symmetrizes the support of the matrices 
Bx,B 2 around the diagonal and saves labour (and terms) in the perturbative expan- 
sions. For the coefficients {p k ,cr k } we make the ansatz 

E P k z k = (63) 



fc=-(*2-l) 

h-l 



z 



E **** = -a-z) h (64) 



*=-i z 



which is in agreement with the support conditions (^8|)-(p0|) if we set l\ = Li, 1% = L 2 
so that the singularity is maximal. Moreover we set 

ci = 1 (65) 

in (|iE|). The normalizations (|63|)-(|6"5|) can be changed; these renormalizations will be 
discussed at the end of this section. 

Inserting (|63|)-(|65D into the Dyson-Schwinger equations (|46|) gives critical coupling 
constants 

A" - i M (66) 
Sf - 9? M (67) 

which can be presented in the explicit form 



/ , , \m/ I m — k 



+ E 

(l>m>m'>fc 



'h\fh(m-l)\ n 2 (m'-l) 
,m) \ m — m' J \ m' — k 



± ... (68) 

The Dyson-Schwinger equations can then be given as two systems of equations for 
the {vk(x)} and {u_ k (x)} 

v k {x) = S k {[u- m (x)),{gi h ' h) }), < fe < Zi - 1 (69) 
u_ k (x) = R- k ([v m (x)],{g^}), 0<k<l 2 -l. (70) 

We first observe that there exist no equations for v_i(x) and Ui(x), so that all 
{v k (x), U- k (x)} are expressed in terms of these functions and their derivatives (plus 
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the {wk{x)}). Second we look at the linear approximation to fl6"9|), ( |7D] 

h-i 



v k lx 



U- 



m=—l 



-'km 



h-1 



E 

k=-l 



and find that its corank is l 2 — 2 (remember l 2 < h) 

Q(ll,h) i s 



= (71) 
= (72) 
The general form of the matrix 



km 



\k+m 



E (n-lX-ltt/f'W 



n=fc+m+2 



' i 2 (« - 2) N 
v n — (A; + m + 2) 



(73) 



Obviously &i m ^ depends on only k + m. Due to the nonvanishing corank (Z 2 > 2) 
we need Z 2 — 2 additional functions which in the perturbative expansion of the 
Schwinger-Dyson equations appear each at order a~ kl . Explicitly this expansion looks 

as 



v 

Vi 

v 2 



\ 



( 



S (0) 



Ml 



-7^(1) 



m?i + a 



-Try ^(2) 



V ^1-1 / 



U-, 



MT 



Mlf-1 

^1 
\ W 2 J 



( 



tflMi 
\ Wlf-1 



+ 0(a" 4 ^ 



(74) 



and correspondingly 



U-2 



Ui 



+ . . 



(75) 



V u -h+i J 

The matrices S^ n ' , have rational numbers as entries. 

On the other hand the Dyson- Schwinger equations imply that these new basis 
functions satisfy (to leading order) nonlinear differential equations 



£> s ({w fe };Mi,u_i) = 0, 



1< s < l 2 - 2. 



(76) 



The introduction of the {wk} is not unambiguous, e.g. one can add to w 2 any linear 
combination of the functions appearing in the same column in ( [74]) as w 2 itself or any 
higher order function such as a" 7 .w'{. The set of all possible choices can be reduced 
by the requirement of gauge invariance (section 6) . 
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Since we have recursion relations for s m , s 



0), < 

ri(m) 
s_i(m + 1) 



the First Lemma implies that all observables may depend only on 

ri(n — ~)s_i(n + |) = 1 + a~ 27 ( , ui(x) + u_i(x)) + a~ 4ry ui(x)v _i(x) 
and other gauge invariant functions. Thus we are led to define as in (||) 

= + + a _27 tti(x)t>_ 1 (x). 

We emphasize that (|7|), ( |50| ) depend on our normalization (|63|) . 

pi = CT-l = 1 



(77) 
(78) 

(79) 
(80) 
(81) 



and must be changed if we renormalize these quantities (see below). Thus we conclude 
that also the functions {wk}i~ 2 must be introduced in such a way that the free energy, 
all expectation values, the leading orders of the differential operators in (|j) and the 
equations ( [T6f) depend only on u and {wk}- 
Now we rescale N such that 



Aa~ 



(82) 



It is obvious from flSTD, flS2]) and (|74l) , ([75]), (|80| ) that to all functions Uk,Vk,u,Wk and 



their derivatives can be ascribed a dimension (or "degree" , "grade" ) under (|8~2"D 

dim{w fc , w} = 2 (all k) (83) 
dim{w fc } = A; + 2 (84) 



dim I — 

ax 



15) 



Let us consider cj as a free parameter in the Dyson- Schwinger equations but maintain 
the normalizations (|63|) , fl64]) . Then the critical coupling constants ( |68|) are replaced 
by 

9k c i9k l° b J 

but the perturbative solution of the Dyson-Schwinger equations is otherwise un- 
changed. The asymptotic expressions ([!]), (@) for -61,-62 remain the same and the 
commutator equation is now 



j2a- n " / p {n \j2 a ~ my Q {m) 



1 (case A) 
a~ 7 (case B) 



17) 
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Using the concept of dimension in the above mentioned sense and the dimensional 
homogeneity of the differential equations ([76]) , (|S?|), we can show that c\ can be tuned 
to one by an appropriate rescaling. 

We can do more. Renormalizing (|63|), (|64|) by replacing the r.h.s. by 



/ 1\' 2 1 
d x z ( 1 J resp. d 2 -(l - z) h 



and keeping c\ as in (|65|), we obtain a new gauge invariant function for which (Q), 
(|80|) is a special case 



(x) = d 2 u 1 (x) + div^i(x) + a 2 ° 'ui(x)v -i(x) . (89) 



This eliminates one degree of freedom in the renormalization. The remaining degree 
of freedom can be eliminated by a rescaling just as c\. Thus we have shown that 
no relevant free parameters are left over in the characterization of the universality 
classes. 

In case (A) the commutator 

[P,Q} = E l s 2 = 2 {r s ,d s }, d=± 

(P = PW,Q = QM) (90) 

does not vanish identically and can therefore be set equal to 1. Now the second 
Lemma asserts that 

r s = 0, se{l,2,... ) / 2 -2} (91) 

whereas 

2r = 1 (92) 
holds by fiat. The equations fl9~T| ) are integrable to 

Jr s dx = h s , s 6 {l,2,...,Z 2 -2}. (93) 

The "constants of integration" h s play the role of external scaling fields. Equation 
( |9"2"| ) can be integrated to 

2 r dx = x. (94) 



The differential equations determine u and all Wk as functions of x and h s . This is 
the standard approach Q. 

On the other hand we have in all cases studied that 

D s (w k ,u)= Jr s dx = 0. (95) 

Since we have tuned only one coupling constant g, the integration constants are 
automatically zero. Of course we can also introduce effective actions, e.g. 

S cS = Tr { Q— + £ ts Q— \ (96) 
12 



for case (A), where 

to ~ x, t s ~ h s (s > 1) (97) 

Thus the known picture (as presented in i, say) is completely reproduced by our 
perturbative approach. 

4. Examples from the series (A) 

In this and the subsequent section we study the differential operator limits of Bi 
and B 2 . For this purpose we need to know the numbers 

P n = ^k n p k (98) 

k 

£ n = £>V* (99) 

n 

U n = Y,k n u k (100) 

k 

V n = ^k n v k (101) 
k 



and the functions 



From (B3D,(B3D results 



0, n < l 2 (102) 

p h = y. (103) 

Pi 2+ i = (1-|)(Z 2 + 1)! (104) 



and 



E n = 0, n<h (105) 

S, a = (-)H! (106) 

= H' 1+1 (l-|)(/i + l)! (107) 
In terms of these quantities the asymptotic expansions of B\ and B 2 are 

CO CO 

a hi Bi = J- q r ( x )d r = a~ n7 Q (n) (108) 



where 



r=0 n=0 



a l ^B 2 = Ps(x)d s = a- n ^P {n) (109) 

s=0 n=0 

q r (x) = \[a-^P r + £ a-C^-W? EMM ] (110 ) 
r. n=0 a n. 
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p.(x) = l[a-^ S s + £ a"^* (111) 

For the series (A) only Q = and P = P^ are of interest. 

In the case (hjh) = (4,3) which describes the critical Ising model, we have only 
one function w — W\ which we define by (exact to all orders) 

uo = u i + V-i + a~ J wi (112) 

It turns out that this definition of w\ can be chosen in all models (see, however, 
([191 ). We obtain from fl7|), (^§ 

C/- = 3wa~^ + 0(a- 2 ^) (113) 

C/j = 3« + 0(a~ 7 ) (114) 

f/ 2 = 0(1) (115) 

V = (§«" + 2w 2 ) a" 27 + 0(a~ 37 ) (116) 

Vi = Awa~^ + 0(a- 2 ^) (117) 

V 2 = 8w + 0(a" 7 ) (118) 

^3 = 0(1) (119) 



This implies 



Q = 8 3 + ^{u,8} + 3w (120) 
P = d A + 2{u, 2 } + 2{w, d} - \u" + 2m 2 (121) 



and 



Di(w,u) = w" + 6wu 

= J ndx = ht = (122) 

These results are equivalent with those known §. 

The case (Zi,^) = (5,4) corresponds to the tricritical Ising model. We have two 
functions wi, w 2 which we define by ( |112| ) and 

u_ x = -5ui + im + 2wi a~ 7 + (6w 2 + v\ + w 2 )«~ 27 (123) 
it follows then from (ff^),([f5|) 

U = (u" + 6u 2 + 4w 2 ) a- 2 ~< + 0(a~ 37 ) (124) 

Ut = 4wi a~ 7 + 0(cT 27 ) (125) 

C/ 2 = 8M + 0(a~ 7 ) (126) 

C/g = 0(1) (127) 

V = (-§ «// - 5uwi) a" 37 + 0(a~ 47 ) (128) 

\4 = (-§«" - 10m 2 - 5w 2 ) a~ 27 + 0(a~ 37 ) (129) 

y 2 = -10w 1 a~ 7 + O(a- 27 ) (130) 

V 3 = -30m + 0(cT 7 ) (131) 

F 4 = 0(1) (132) 
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The two differential operators come out as 

Q = d 4 + {2u, d 2 } + {2w u d} + Aw 2 + 6m 2 (133) 



2 1 

After the substitutions 



P = 8 5 + -{u, d 3 } + -{<&!, d 2 } + -{w 2 - \u" + 2u 2 , 0} + 5u Wl (134) 



1 . . 

u — > —-u (135) 

W\ — > -w (136) 

Aw 2 + 2u 2 -> v (137) 

these differential operators and the differential equations (|92|), (|95|) can be shown to 
be identical with those of Ginsparg et al. 0. 

Now we come to the case (h,l 2 ) = (6,4) which was not studied before. In this 
case there is a common divisor of l\ and l 2 : 2. Once again we can define w\ and w 2 
as in (|112|) , (|123|) but we prefer to define w 2 h y 



U-i = —5u\ + V-i + 2w\ a 7 + (|w 2 — \v % _ x — ?>u\V-\ + w 2 )a 27 (138) 

because this minimizes the terms in the differential equations that we have to give 
explicitly. The differential operators are 

Q = (d 2 + 2u) 2 + 2{ Wl ,d} -Aw 2 (139) 

P = (d 2 + 2u) 3 + 3{w 1 ,d 3 }-3{w 2 ,d 2 } + {-lw'( + 6uw ll d}-2w 2 +6w 2 1 -12uw 2 (140) 
By explicit calculation we found that 

P = Ql (141) 

Moreover we get the differential equations 

h 2 = D 2 (w 1 ,w 2} u) 

= —w[^ — 8uw" — 12u'w[ — Au'wi — 2Aw\w 2 (142) 
hx = D 1 (w 1 ,w 2 ,u) 

= w { 2 ] + 8uw' 2 ' + Au'w' 2 + I2w 2 2 - lOwiw" - 5w[ 2 - 2Auw 2 (143) 

x = y*[2r + \r 2 + 2u'(D 2 (w 1 , w 2 ,u) - h 2 )]dx 

= SwtwZ - 8w[w' 2 - 16wf - 2h 2 u (144) 

To make 2r integrable weJiad to add the term 2u'(D 2 — h 2 ). But this is also necessary 
in the (5,4) example (see 0). 

5. Examples from the series (B) 

The simplest non trivial example is the symmetric case 

(M 2 ) = (3,3) 
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There exists one function W\ — w which we introduce by (|112| ). In this example we 
obtain from perturbation theory 



- p(°) =Q(0) = dP + ±{ U ,d} + Zw 



and 



(145) 
(146) 
(147) 



Di(w, u) = w" + 6wu = 

PW and are separately not gauge invariant. 

Thus we have reproduced the differential operators and equations ( |120| )- (|122|) 
The commutator 



Y,a~ nj P {n \Y, a ~ m ^Q (m) = [~ p{1) -Q {1 \Q {0) ] a" 7 + 0(a" 27 ) 

n m J 



(148) 



of a" 7 in [P 2 ,-Bi] to 

(149) 



gives an additionnal factor a 7 that enlarges the total power 

3+3+1=7=4+3 

as is necessary for the universality class [4,3]. 
Now we turn to the rather complicated case 

(M 2 ) = (M) 

Again we introduce w\ = w by ( |112| ). The differential operators are 
Q (0) = 8 3 + ^{u,d} + 3w 

Q« = -^-^{u 1 + 7v^,df}-{w,d} + ^(u , (+13v'U)+3u 2 1 -u 
p(°) = (QW)2 

p« = 2d 7 + ^{l7u 1 + llv- ll d 5 } + 7{w,d A } 

+-{-i(281< + 161^0 + 43m? + 56«iv_i + 19v\, d 2 } 
' ' "" , "" ' 47^)-|(155ui< + 119n>_i + 107rf 1 + 71w_i^ 1 

/ / , rrf.J \2\ 



(150) 

(151) 
(152) 



+{i(83<' + ... N , 

-!(15K) 2 + 22^ 1 + 7(0 



+2Qu\ + 42m 2 v_i + 24«it; 2 _ 1 + 8vl t + 3fV , d} 

+—w"" + 3w" - -w'(49u[ + 31^) - 7^(15< + 3v'L lt 
6 2 4 

+18u>(3ti 2 + 4miw_i + 2w 2 x ) 



(153) 
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Both P« and are gauge dependent. We consider the commutator 



J2 a -n 7p (n) ^ J2 a~ m7 Q (m) 

n m 

n m m 

[ p(D _ Q(o)gd) _ q(Dq(o) } q(o) ] a -7 + 0(a- 27 ) 

The first differential operator in ( |154| ) can be calculated and is 

p(i)_ ( Q(0)Q(i) + g(i)g(0) )=3( Q(0) ) I 



(154) 



(155) 



The factor 3 can be eliminated by rescaling as was discussed in section 3. The 
additional factor a -7 changes the total power of a -7 in [B 2 , Bi] to 



6+3+1=7+3 



(156) 



as is necessary for the universality class [7, 3] . Thus this class has been proven to 
arise. 

We remember that one-matrix models lead to universality classes [2m + 1,2], and 
we could imagine that such classes can also be produced from two-matrix models with 
singularities 

(/ 1 ,/ 2 ) = (2m,2) 

However , if we insert an ansatz for such singularity into the Dyson-Schwinger equa- 
tions, the perturbative solution cannot be made to satisfy the commutator equation 



[B 2 ,B 1 



1 



(157) 



at any order in a~ 7 as we shall show in the sequel. 

Of course a quadratic potential V 2 can be eliminated from the partition function 
by performing a Gaussian integral. In the Dyson-Schwinger equations this comes out 
as follows. If 



A 1 + cB 2 = V&BJ 

A. + cB, = gf ) +gf ) B 2l (gf^O) 

then B 2 can be eliminated by substituting (|159|) into ( |158| ) 

( ' fi {g { ? - cflO + V[{B X ) 



(158) 
(159) 



A, + —A 2 

92 



(2) 
92 



= V(B 1 ) (160) 

Now Ai is strictly lower and A 2 strictly upper triangular and under transposition 

£?! -> B\ (161) 
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they exchange their role in ( |160| ). In the one-matrix model there is only one orthogonal 

system {IT m (A)} (section 2) implying that B x is symmetric. In this case ( |160| ) reduces 
to 

A x +Aj — V'{B X ) (162) 

Via the antisymmetric matrix 

C = A X - -V'iBt) (163) 

and the double scaling limit 

B x -> a~ 27 Q (164) 

C -> a~ (2m+1)7 P (165) 

the universality class [2m + 1,2] is obtained. 

Instead we solve the Dyson-Schwinger equations ( p.58| ), (|159[) perturbatively as 
they stand. First we set 

L x = l x = 4 , m = 2 (166) 
Our algorithm simplifies by the equalities to all orders 

u Q (x) = v Q (x) (167) 
u-i(x) = v-i(x) (168) 

Due to the fact that no functions {wk} are needed, we obtain only even powers in the 
expansions ([!]), @ 

oo 

B i -» a ^ £ a- 2n7 Q (2n) (169) 

n=0 

oo 

B 2 -> a" 47 ^ a -2n 7 p(2n) ^yg) 



We find 



n=0 



Q(0) = 9 2 + 2w (171) 

Q (2) = ^ + \{u,d*}-\u"-\{u x -v„ x T (172) 

p(0) = ( g(0 }) 2 (173) 

P(2) = ^ 6 + 7^{ 4lM i + 11 ^-i^ 4 } 

+{- §(23< + 3^0 + 16u 2 + lhu x v- x + <9 2 } 

+i(5< + ^ 1 )-l(15< + 9^ 1 ) 

+ 23u x u" + 5^1^^ — Sf-iV/ + 7f_if" 1 ) 
8 

+ ^(27m? + 67w 2 f_! + 37m!!; 2 ! - 3v 3 _ 1 ) (174) 



An almost trivial observation is that all P^ 2n \ are symmetric operators. 

Consequently the commutator 



£ a -2n 7 p(2n) ^ £ a -2 m7 g(2 m ) 
n m 
[ p(2) _ Q(0)Q(2) _ g(2)g(0) ) g(0) ] fl -2 7 + 0(fl -4 7) = 

E{ri 2) ,a fc }a^ + 0(a- 4 ^) (175) 

k=0 

is antisymmetric 

r f = r f = r f = (176) 
This symmetry argument is valid to any order (and for any m) 

rf n) = , all n (177) 

and therefore the commutator relation cannot be satisfied. 

Only as a side remark we add that the three differential equations 

= rf ] = rf ] = (178) 

ux = u_! (179) 
Bi = Bj (180) 



have the unique solution 
which implies 



6. Expectation values and gauge invariance 



The basis functions {wk} have hitherto been introduced (see ( |112| ), ( |123| ),( |T3S| )) by 



the requirement that the relevant objects, namely differential operators and equations 
are gauge invariant. This fixes them only at low orders. Expectation values are a 
simple device to extend this definition to higher orders. 

For any two- matrix model we obtain for a finite iV using (|31|), fl37D,([48D and (|77|), 
(|78|) that the partition function can be expressed as 

N-l 

Z = C(N)N\(s ~s r N n (n(N — k — l)s_ x (JV - k)) k (181) 

k=i 

so that the free energy is 

F = \ogZ = log(C(N)N\)-Nlog(s s ) 

+ a" 27 (l - a 2 z)~ 2 [ Z dx(z - x) log(l + a'^uix)) 

+ corrections from Euler's summation formula. (182) 
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The latter corrections in (11821) form an infinite series. The function u(x) increases for 



x — > oo and one can convince oneself that 

a- 2 ^u{a- 2 ) = 0(1) (183) 



so that the perturbative expansion of the integral in ( |182| ) fails at the lower end. The 
corrections due to Euler's summation formula contain terms 

u (n \a- 2 ), n>0. (184) 

Therefore we set 

F(z) = F(0) + (F(z) - F(0)) (185) 

and expand only the difference perturbatively. A similar caveat must be applied to 
the expectation values. 
Expectation values 

(G/(M«,M^)) (186) 

with respect to the partition function (|25| ) or ( |31~D can apparently be obtained by 
parametric differentiation, if G is a polynomial of the variables 

Tr(M«) a , Tr(M^) b , Tr(M^M^). (187) 

After reduction of the partition function by means of Mehta's theorem, these variables 
correspond to 

£^> E A ^> (A i = A? ) ,A*i = \ p,) ). (188) 

i i i 

First we show that the expectation value of the third variable with a reduced partition 
function is not gauge invariant. 

Let us denote the trace over the upper left N x N submatrix of an infinite (IN x IN) 
matrix by Tr^- Then insertion of 

i 

into (|31~D gives for fixed iV 

Ti N {B l B 2 ) 



which can be evaluated to yield 
N 



% + h 



(/l + /2 - 1)! (/ 1 / 2 -/ 1 -^ 2 ) + l 



a~ 7 



liH 2 \ 

f_ 2 dx | Y,[p-kVk(x) + a k u„ k (x)} \ + 0(a- 37 ) (189) 



for any maximal singularity I2). For (Jx, l 2 ) = (4, 3) the integral is to leading order 
(using still ([1121) ) 

Uu x {x) - 43v_i(x) - 17wi(x)a -7 (190) 
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which is not gauge invariant. Therefore only the first two expressions in ( |188| ) make 
sense. _ 
Expectation values of monomials can be reduced easily (see Ej) 

(Tr(MW)°Tr(M«) 6 ) = Tr N (B«)Tr N (Bl) +Tr N (BfII N B b 1 ) (191) 
(Tr(M (1) ) a Tr(M (2) ) b ) = Tr N (B*)Tr N (B b 2 ) + Tr N (B°Tl N B b 2 ) (192) 

where II at is the projection matrix 

, x / S rs if r > N and s > N , , 

(^N)rs ~ Q elge _ (193) 



Of course for the simple cases of one variable ( 187|) we have from ( |191|) , (|192|) 



(Tr(M (1) )°) = Tr 7V (5 1 ) a (194) 
(Tr(M^) b ) = Tr N (B 2 )». (195) 

Consider (|194| ) with a = 1. Then 

Tr^Bi = Np — a -7 J ^ dxuo(x) 

+ corrections from Euler's summation formula. (196) 
We recognize that instead of (|112|) we better use 



u (x) = u\ + v _i + a 7 Wi + a 27 uiV-i 

= u + aT'wx. (197) 

The same argument can be applied to 

Tr^B 2 = Nan — a' 1 \ dxv (x) 

Ja- 2 

+ corrections from Euler's summation formula. (198) 

so that vq must depend only on u and w\. Indeed we find for the maximal singularity 
(4, 3) from ([197]) 

v = u- Wl a-~< - + u 2 )a- 2 ~ f + + u Wl )a- 3 ' ( + 0(a~ 47 ) (199) 

and for (3, 3) 

Vq = u — - (±u" + u 2 )a~ 2 ^ + 2{\w'{ + u Wl )a-^ + 0(a" 47 ). (200) 



Finally we consider the cases a = 2 in ( |194| ) and b = 2 in (|195| ) 
To sum up the results we introduce abbreviations 



<j>x(x) = u{x) (201) 

falx) = wx(x) (202) 

fa(x) = ±u"(x) + u 2 (x) (203) 

(f) 3 (x) = \w'{{x) + u(x)w 1 (x) (204) 
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$ n (z) = f <f) n (x)dx. (205) 

Ja~ 2 

We write for (PI , QT95D 

7iiV + 7o + E c « $ nW«" n7 (206) 

n>l 

and obtain the coefficients in the following table 



sing 


a 


b 


7i 


7o 


Cl 


c 2 


c 3 


c 4 


(4,3) 


f 




-3 





-1 


-1 










2 




15 


-3 


4 


2 


-1 


22 
5 






1 


-4 





-1 


1 


1 


2 
5 






2 


28 


-6 


6 


-4 


-5 


-10 


(3,3) 


1 




-3 





-1 


-1 










2 




15 


-3 


4 


2 


-1 


2 






1 


-3 





-1 


1 


1 


-2 






2 


15 


-3 


4 


-2 


-3 
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Table 1. Expansion coefficients for the ex- 
pectation values of ( |194j ), ( |195| ) for the sin- 
gularities (4,3) and (3,3). 



The general expression for 71 is 

(-l) a ( ah ), respectively 
7i= V 1 (207) 

1 (-l) 6 ( 6 f). 
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